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A comparative study of overlap and staggered fermions in the Schwinger 
model* 
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We investigate the validity of the square rooting procedure of the staggered determinant in the context of the 
Schwinger model. We find some evidence that at fixed physical quark mass the square root of the staggered 
determinant becomes proportional to the overlap determinant in the continuum limit. We also find that at fixed 
lattice spacing moderate smearing dramatically improves the chiral behavior of staggered fermions. 



1. Introduction 

Recently, unquenched lattice QCD calculations 
with improved staggered fermions have had re- 
markable success in reproducing a variety of phe- 
nomenologically interesting quantities Pj. In 
these calculations, the determinant of a single 
fermion is obtained as the fourth root of the stag- 
gered determinant. It is a priori unclear whether 
there exists a local operator D describing a single 
fermion flavor with 



det (£)) oc det (-Dstag) 



1/4 



(1) 



and it has been shown 2 that a naive guess leads 
to a nonlocal operator. Furthermore, the lack of 
an (exact) index theorem for staggered fermions 
raises the question how far into the chiral regime 
one can push calculations at fixed lattice spacing. 

Given this situation, we decided to study the 
behavior of staggered fermions in the simple and 
well known Schwinger model (QED2). We focus 
on observables that - like the determinant - are 
obtainable from the Dirac spectrum alone. As a 
point of reference, we compare the staggered data 
to results obtained with overlap fermions, which 
are known to be free of conceptual problems. We 
also investigate UV-filtered staggered and over- 
lap fermions, which are technically realized by 
APE-smearing the gauge backgrounds. For the 
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U(l) gauge group smearing consists of taking a 
weighted average of the phases of the original 
link and the staple. Due to the twofold fermion 
doubling in 2D, the Nf = 1 case is obtained 
by weighting the gauge configurations with the 
square root of the staggered determinant. Tech- 
nical details of our simulations can be found in 
0, related studies in QCD are reported in [4l5j . 

2. Chiral condensate 

Due to the Mermin- Wagner theorem, the chiral 
condensate of the continuum massless Schwinger 
model vanishes for Nf > 2 For Nf = 1 the 
formation of a condensate does not imply sponta- 
neous symmetry breaking. In the massless case, 
its infinite volume value is analytically known [7] 
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Away from the chiral limit one expects an ad- 
ditional contribution to the condensate propor- 
tional to the quark mass and at most loga- 
rithmically divergent in the cutoff. Also, due 
to the exact /remnant chiral symmetry of over- 
lap/staggered fermions, the chiral condensate 
renormalizes multiplicatively and due to the di- 
mensionful coupling [e] = [a]~^ all renormaliza- 
tion factors are Z = 1 + 0{a^e^). 

In fig. n we plot the Nf — 1 chiral condensate 
versus the quark mass, as obtained from stag- 
gered and overlap fermions. Our results for the 
unfiltered operators (dotted lines) indicate that 
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Figure 1 . {ipi^) /e of staggered and overlap 
fermions with Nf = 1 on a 242 lattice at /? = 7.2. 
The star denotes the analytically known value of 
the continuum chiral condensate from |2|. 
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Figure 2. Comparison of staggered versus over- 
lap spectra for three typical configurations (3 left 
columns) and one configuration with undecided 
topology (right column). 



overlap fermions show a nice chiral behavior while 
staggered fermions behave qualitatively wrong in 
the chiral limit. By contrast, the UV- filtered 
staggered fermions (full line) exhibit a much bet- 
ter behavior down to comparatively small quark 
masses. They still show, ultimately, a qualita- 
tively wrong behavior in the chiral limit, but the 
mass at which it sets in is drastically reduced. 

3. Spectral mimicry 

A finite condensate in the chiral limit is related 
to the existence of exact zero modes. The dra- 
matically improved behavior of {'tfjijj)/e for UV- 
filtered staggered fermions down to fairly small 
quark masses therefore suggests that they develop 
near-zero modes on topologically nontrivial gauge 
configurations and that very small quark masses 
are needed to expose these modes as not being 
true zero modes (see [H] for a more comprehensive 
discussion of this point). It is thus interesting to 
directly compare staggered and overlap spectra 
on individual configurations. 

Fig. 13 shows a comparison of the IR part of 
massless staggered and overlap^ spectra on some 
selected configurations on a 20^ lattice at /? = 4. 

-"^For the overlap operator, we plot the (purely imaginary) 
eigenvalues A = (1/A — l/2p)~^ of the chirally rotated 
operator (see |21 for more details). 



The two plots on the left show typical configu- 
rations with topological charge q = and \q\ = 1, 
respectively. The unfiltered staggered spectra do 
not resemble the overlap spectra and it is difficult 
to see the qualitative difference between q = 
and |g| = 1. After moderate UV- filtering how- 
ever, the eigenvalues of the staggered operator 
form near-degenerate pairs which sit close to a 
single overlap eigenvalue. In particular, in the 
l^l = 1 case a pair of eigenmodes moves very close 
to the real axis, confirming the expectation about 
the existence of near-zero modes that we had from 
the chiral condensate analysis. 

The third column of fig. [5] shows a typical con- 
figuration with higher topological charge (jgl =4). 
Here, after 3 smearing steps only 3 pairs of eigen- 
modes have managed to come close to the real 
axis and the fourth one is still further out. 

The last column of fig. 12 presents a selected 
"worst case" configuration on which the topo- 
logical charge varies repeatedly under subsequent 
smearing steps. Qualitative resemblance between 
the two types of spectra is vague at best and there 
is no evidence for a pair of staggered eigenvalues 
moving close to the real axis. Such configurations 
are rare already at /3 = 4.0. 



3 



4. Determinant ratio 

The qualitative similarity between the IR part 
of staggered and overlap spectra on individual 
configurations leads to an interesting suggestion 
regarding the local operator D in If eigen- 
value pairs become truly degenerate in the con- 
tinuum limit and topologically undecided config- 
urations arc arc strictly O(a^) effects, it could be 
that I?ovcriap is a local operator with 
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ocdet(i:iovcriap) + 0(a^) . (3) 



Fig-Elshows a scatter plot of log(det(Dstag)^^^) 
vs. log(det(Dovoriap)) for different lattice spacings 
and UV-filtering levels at one fixed physical quark 
mass. One can clearly see that there is a tendency 
for the determinants to be proportional in the 
continuum limit and that the proportionality is 
better developed for the UV-filtered operators. 

With a careful scaling study one might be able 
to test Q. The real goal is, of course, to see 
whether such a relation holds in QCD. Prelimi- 
nary studies in this direction can be found in 
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Figure 3. Log-Log scatter plot of the over- 
lap determinant versus the square root of the 
staggered determinant at 4 lattice spacings and 
3 UV-filtering levels and fixed fermion mass 
m/e = 0.035. Every plot shows 1000 decorrelated 
quenched configurations at fixed physical volume. 
The individual plots are offset for better visibility. 
In all cases the slope is consistent with 1. 



